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A topographical study on an ensemble of height profiles obtained from atomic force microscopy techniques on various 
independently grown samples of tungsten oxide WO3 is presented by using ideas from percolation theory. We find 
that a continuous 'geometrical' phase transition occurs at a certain critical level-height 5 C below which an infinite 
island appears. By using the finite-size scaling analysis of three independent percolation observables i.e., percolation 
probability, percolation strength and the mean island-size, we compute some critical exponents which characterize the 
transition. Our results are compatible with those of long-range correlated percolation. This method can be generalized 
to a topographical classification of rough surface models. 
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The growth of rough surfaces and interfaces with many of 
their scaling and universality properties has attracted the at- 
tention of statistical physicists over the last three decadesA. 
One of the less studied subjects is the one concerned with the 
topographical properties of the self-affine surfaces. In this let- 
ter, we present a systematic investigation of percolation tran- 
sition on an ensemble of experimentally grown WO3 surfaces. 
This gives some 'geometrical' characteristic exponents which 
well describe a geometrical phase transition at a certain level- 
height below which an infinite cluster-height (or island) ap- 
pears. The results of our work are closely related to the statis- 
tical properties of the size distribution of the islands and their 
coastlines together with their fractal properties. 

The tungsten oxide WO3 surface is one of the most interest- 
ing metal oxides. It has been investigated extensively because 
for its distinctive applications such as electrochromio^— , 
photochromies-^, gas sensor— ~—, photo-catalyst^, and photo- 
luminescence properties-^. Many properties of the WO3 thin 
films are related to their surface structure as well as to their 
surface topography and statistics such as grain size and height 
distribution. These properties can also be affected during the 
growth process by either deposition method or imposing ex- 
ternal parameters on the grown surfaces such as annealing 
temperature-^. 

In order to have an ensemble of WO3 height profiles, 34 
samples were independently, and under the same conditions, 
deposited on glass microscope slides with an area 2.5 cm x 2.5 
cm using thermal evaporation method. The deposition system 
was evacuated to a base pressure of ~ 4x 10~ 3 Pa. The thick- 
ness of the deposited films was chosen to be about 200 nm, 
and measured by stylus and optical techniques. Using atomic 
force microscopy (AFM) techniques, we have obtained 300 
height profiles from the grown rough surfaces with resolution 
of 1/256 /im in the scale of 1 /imx 1 /jm (an example is shown 
in Fig. [TJ. The AFM scans were performed in various non- 
overlapping domains (10 images from each sample) from the 
centric region of the deposited samples in order for the profiles 
to be statistically independent. 



The percolation problem-! 15 j s an example of the simplest 
pure geometrical phase transitions with nontrivial critical be- 
havior, and it is closely related to the surface topography. Let 
us suppose a sample of height profile {h(r)} is a topograph- 
ical landscape. Now, let us imagine flooding this landscape 
and coloring the parts above the water level white and the rest 
black. If the water level is high, there will be small discon- 
nected islands, while if it is low, there will be disconnected 
lakes. There is however a critical value of the sea level for 
which there is one large supercontinent and one large ocean. 
As long as the original height profile has a gaussian distribu- 
tion with only short-range correlations, it is believed that the 
large-scale properties of the coastlines correspond to standard 
percolation cluster boundaries and thus should be described 
by the theory of Schramm-Loewner evolution (SLE)!£ (seeii 
for a review of SLE). However, the height profiles of WO3 
surface are quite different due to the relevant contribution of 
long-range correlations-^. Hence, the corresponding coast- 
lines of the height profiles belong statistically to a different 
universality classic. 

To define islands on a WO3 surface, a cut is made at a 
certain height h s = (h(r)) + <V(IM r ) - (M r ))] 2 > : = 0, 
where the symbol (••) denotes spatial averaging. Each island 
(cluster-height) is then defined as a set of nearest-neighbor 
connected sites of positive height. We show that there is a 
critical level height denoted by the dimensionless parameter 
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FIG. 1 . (Color online) AFM image of WO3 thin film in scale 1 /zm 
x 1 fim with resolution of 1/256 /im. 
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FIG. 2. (Color online) Main: probability P s for the presence of a 
spanning island as a function of 5 measured for different lattice sizes 
L. The curves cross at a critical level height S c = —0.20(1). Inset: 
data collapse for the P a curves of different L with v — 1.90 and 
S c = -0.20. 



S = 5 C , at which a continuous percolation transition occurs. 

To study the finite-size scaling (FSS) properties and mea- 
sure percolation observables which characterize the critical 
behavior, we consider a box of dynamic size L x L from the 
centric region of each original AFM sample. For each lattice 
size 16 < L < 256, each percolation observable was aver- 
aged over all number of AFM samples. 
The first quantity we measure is the probability P s that at each 
level height 6, an infinite island spans two opposite boundaries 
of the box in just a specific direction, say y-direction. Ideally, 
the curves obtained for different lattice sizes cross at a single 
point, marking the critical level height S c , As shown in Fig. [21 
the measured curves cross at S = 5 C , implying that the scaling 
dimension of the percolation probability P s is zero. 
According to scaling theory 2 ^, one expects that all the mea- 
sured curves should obey the scaling form 

P s (6) = P s [(6 - 8 c )lV»], (1) 

where the exponent v characterizes the divergence of the cor- 
relation length £ (proportional to the spatial extent of the is- 
lands) near the percolation threshold, £ ~ \S — S c \ ~ v . 
We measure the values of the exponent v and the crossing 
point of the curves S c by utilizing the data collapse. The qual- 
ity of the collapse of the curves is measured by defining a 
function S(is,S) of the chosen values v and S (the smaller S 
is indicative of a better quality of the collapse — se o 21,22 and 
the appendix o£l). We find its minimum S m i n ~ 1.78 for 
v = 1.90(30) and S c = -0.20(1). Inset of Fig. [2] illustrates 
the collapse of all the P s curves, within the achieved accuracy, 
onto a universal function by using the estimated values for v 
and 5 C . 

Percolation strength P^, which measures the probability 
that a point on a level height 5 belongs to the largest island 
(or equivalently, the fraction of sites in the largest cluster), is 



another quantity defined as order parameter in the percolation. 
As shown in Fig. [3] we have computed P^ as a function of S. 
We find that the data follows the scaling form 

P 00 (6) = L-PP 00 {(S-5 C )L 1 / V ], (2) 

Our best estimate for the exponent (3 ((3 = fijv in percolation 
theory) is /3 = 0.12(4) which was obtained by optimizing the 
quality of the collapsed data. The data collapse is shown in 
the inset of Fig. [3] 

Another independent observable is the mean island size 

where n s denotes the average number density of islands of 
size s, and the prime on the sums indicates the exclusion of 
the largest island in each measurement. 
As presented in Fig. |U the obtained curves x($) f° r differ- 
ent lattice sizes have their maximum around the critical level 
5 C . We find that x(S c ) ~ ^ with 7 = 1.70(3) (7 = 7/V 
in percolation theory)— see inset of Fig. SJa). By using the 
exponents 7 and v, it is possible to achieve a data collapse 
according to the scaling form 

X {5)=L\[{5-5 c )L l /»], (4) 

which is shown in Fig. Sfb). 

A point which remains is determination of the universal- 
ity class which the observed percolation transition belongs to. 
The values we find for the exponents are obviously not those 
of short-range correlated percolation. In order to see whether 
they fit with long-range correlated percolation, we calculate 
the two-point correlation function G(x,x') which is defined 
as the probability that two points of distance r = |x — x'|, 
at the critical level height S c , belong to the same island. The 
best fit to our data shows that G is a decreasing function of 
the distance with a power-law behavior G(r) ~ r~ v and 
r] = 0.34(2). 
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FIG. 3. (Color online) The percolation strength of the islands as 
a function of 8 for different lattice sizes L. Inset: collapse of the data 
in the suitably rescaled scales. 
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FIG. 4. (Color online) (a) mean island size x as a function of S for 
different lattice sizes L. Inset: \ at the critical level cut 5 C — —0.20 
as a function of L. (b) data collapse for \ f° r different lattice sizes, 
according to Eq. (@. 



In the past, several papers have dealt with percolation with 
long-range correlation a 1 8 - 24 ^ in which correlations decrease 
with increasing r. A completely different percolation with 
long-range correlations^ has also been proposed by Sahimi 28 
in which correlations increase as r does22. From the above 
references, our results are in agreement with 26 in which the 
value of v is consistently lower than the prediction v = 2/rj 
given byi£, for r\ < 1 in two dimensions. 
In order to examine whether our exponents satisfy well-known 
scaling and hyperscaling relations, we have also computed the 
fractal dimension of the islands at the critical level S c , and 
found to be df = 1.83(3). It is thus straightforward to see that 
our obtained exponents satisfy, within the statistical errors, the 
following scaling relations for d = 2, 

d f =d-p=±(d + j) = ^(d + 2-ri). (5) 

In conclusion, analysis of independent percolation observ- 
ables on an ensemble of height profiles obtained by AFM 
images of WO3 surfaces, revealed a continuous geometrical 
phase transition at a certain critical level height S c . We com- 
puted some critical exponents which can be regarded as topo- 
graphical characteristics of the WO3 surfaces. This method 
may lead to a topographical classification of self-affine rough 
surfaces by computing the exponents v, j3, 7 and 7/ for differ- 
ent growth models. 
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